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I. INTRODUCTION 



The six-vertex model was first introduced as a model for the residual entropy of ice and for related ferroelectric 
transitions |^ but more recently, several other physical applications are being found. In the body-centered solid-on- 
solid models, the six-vertex configuration is mapped to the surface configuration of the fee and bcc faces and the 
free energy of the six-vertex model as a function of the horizontal and vertical fields depicts the equilibrium crystal 
shapes (ECS) Also, on the stochastic surface where the vertex weights satisfy a certain relation, the transfer 

matrix of the six- vertex model can be regarded as the transition matrix of probabilistic cellular automata describing 
the dynamics of a driven lattice gas system and a (1+1) dimensional surface growth model j4|-|^. The five- vertex 
model, which is a special asymmetric limit of the six-vertex model, can be viewed as models for an interacting domain 
wall system 0, interacting dimers jsj, and certain types of crystal surfaces In these applications, it is necessary 
to consider the asymmetric six- vertex (ASSV) model in which there are non-zero horizontal and vertical fields which 
break the arrow reversal symmetries. It becomes the symmetric six- vertex (SSV) model when the fields are zero. The 
Heisenberg XXZ chain is a closely related problem. The Hamiltonian for the XXZ chain can be obtained from an 
anisotropic limit of the transfer matrix of the six- vertex model |^ . The vertical field corresponds to the magnetic field 
which couples to the z component of spins and the horizontal field corresponds to the asymmetry in the two hopping 
rates. The XXZ chain is called symmetric (asymmetric) if the two hopping rates are same (different). 

The phase diagram and the nature of the phase transition in the ASSV model are well known [|l 10|. The phase 



diagram consists of the ferroelectrically and anti-ferroelectrically ordered phases and the disordered phase. Especially, 
the disordered phase shows interesting scaling properties. It is the critical phase with continuously varying critical 
exponents and is described by the central charge c=l conformal field theory (CFT) in the continuum limit [ pd| . Since 
the c=l theory plays a basic role in the theory of two dimensional critical phenomena, there have been many works on 
the scaling properties of the six-vertex model. However, most of these works are confined to the SSV model |jl^ or the 
symmetric XXZ chain with or without the magnetic field ||l^-p^ and to the special case of the five- vertex model , 
while the critical properties of the general ASSV model are also of interest. One purpose of this work is to fill this 
gap. In this paper, we investigate the critical properties of the ASSV model through the finite-size scahng (FSS) 
studies on the transfer matrix spectra. Our work also covers the critical properties of the asymmetric XXZ chain. 

The FSS is very useful in studying the critical behaviors of two-dimensional statistical systems on lattice [0 . When 
the system is in a critical phase which possesses the conformal symmetry, it imposes strong restrictions on the FSS 
form of the eigenvalue spectra of the transfer matrix. Let T denote the row-to-row transfer matrix for a system with 
N columns and M rows. We assume periodic boundary conditions unless stated otherwise. The energy of the level a 
defined by Ea = — In A^, A^ being the a-th eigenvalue of T, is expected to follow the scaling form 

Ea^Nf+ ^C" {he + ha-{^)- {ha -K)+0 (^1) (1) 

where / is the bulk free energy (in units of /c^T), C (C") is the real (imaginary) part of the complex anisotropy factor 
C, = C,' -\- iC," , c is the central charge, {ha, ha) are the conformal dimensions associated with the level and o{x) stands 
for terms smaller than x; limjj^o o{x)/x = 0. The level whose energy follows the scaling form of Eq. (Ill) is associated 



1 



with a primary operator and its descendants of the corresponding CFT. When the one lattice-unit translation by 
T is incommensurate with the periodicity of the underlying model, an imaginary term of 0(1) may appear in the 
right-hand side of Eq. (|l|) Q. The toroidal partition function (TPF) Z is defined as the 0(1) part of the partition 
function: 



Z= lim (2) 

M/Ar=fixcd " 

where the sum is over all levels. If we use the scaling form of Eq. (|^), it takes the form 

^=(qq)-*Eq'°^'° (3) 

a 

where the nome q = e^'^*'^ with t = r' + ir" = -^C and q is the complex conjugate of q. The complex parameter 
T is the modular ratio of the torus on which the corresponding CFT is defined and specifies how the N x M lattice 
should be deformed to make the system isotropic in the continuum limit. The TPF contains the whole information 
on the spectra or the operator content of the model and enjoys the modular invariance (or more generally, modular 
covariance) properties which together with the conformal invariance principle is often sufficient to determine the form 
of Z The c—1 CFT consists of three isolated points and two one-parameter families ||2^]. The latter two describe 

the critical eight-vertex model or the Ashkin-Teller model and the SSV model or the Gaussian model compactified 
on a circle, respectively. They are related by a duality transformation. The TPF of the c=l theory corresponding to 
the SSV line is the so-called Coulombic partition function given by, for M and N even. 



where r]{q) is the Dedekind eta function, 

oo 

,y(q)=qi/24]^(l_q") (5) 

n=l 

and the conformal dimensions are given by 



A™.n = i (7I - Van 
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(6) 



Here g is the Gaussian coupling constant and the integer indices m and n label the spin-wave excitation and the 
vortex excitation, respectively. The Gaussian coupling constant g is defined in such a way that it takes the value 1/2 
for the free fermion theory and is related to Kn of Ref. [|l| by g = 27ri^/j. However, the CFT being a bootstrap 
theory, it does not tell us how g and r are related to the lattice model parameters. To obtain that information, one 
needs to rely on the FSS analysis. 

The transfer matrix for the six-vertex model or equivalently the Hamiltonian of the XXZ chain is diagonalized 0| 
by the Bethe ansatz method for general boundary conditions. In this paper, we present a method of calculating 
the finite-size corrections of arbitrary low-lying energies in the whole parameter space of the ASSV model. Starting 
from the Bethe ansatz solution for the system of width iV, we derive a systematic expansion in 1/N of the energies 
assuming certain analyticity properties of the phase function which is introduced in Sec. The method is very 
similar to that used in [p^ which considered the symmetric XXZ chain in a magnetic field but is generalized to be 
applicable to the general cases. From the expansion, we find that the whole critical phase of the ASSV model is 
also in the Gaussian- model universality class with c = 1, where the coupling constant g is obtained from the solution 
of an integral equation. And we also find that though the susceptibilities are non- universal, i.e., they depend on 
non-universal parameters, the Hessian of the free energy is simply given by 



H = 



91/ d^f 
dih dhdv 



dvdh 



7rg 



(7) 



and hence depends on the model parameters only through g. Here, h and v are the horizontal and vertical fields, 
respectively, in units of fcsT. Since we find the FSS form for a general class of energy levels, we are able to construct 
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the TPF for the ASSV model with the periodic boundary conditions exphcitly. The resulting expression is given by 
the modified Coulombic partition function 



^= E e— V"."-^q^"."-^ (8) 

where a and f3 as defined in Eq. (^2|) are the two mismatch parameters which account for the incommensurability of the 
lattice with the mean distances between down arrows and left arrows, respectively. The ASSV model under the twisted 
boundary conditions (explained in Sec. ||) is equivalent to the ASSV model under the periodic boundary conditions 
with modified fields. Using this relation the TPF for the ASSV model under the twisted boundary conditions is 
also derived, which confirms the conjectured TPF for the SSV model under the twisted boundary condition in the 
horizontal direction . 

This paper is organized as follows. In Sec. |^ we give a brief review of the transfer-matrix formulation and its Bethe 
ansatz solution of the ASSV model. The classification scheme of the Bethe ansatz solutions is presented. The relation 



between the ASSV model under the periodic and the twisted boundary conditions is also discussed. In Sec. Ill, under 



a certain assumption we derive a summation formula which converts a sum over functions of rapidities for general 



levels into a finite-size expansion. The assumption turns out to be the sufficient condition for the criticality in Sec. IV. 
We obtain an integral equation for the expansion coefficients and find that the finite-size correction terms are related 
to the partial derivatives of the bulk term with respect to the horizontal and vertical fields. Some of the details of 
calculations are relegated to App. In Sec. IV, applying the summation formula to the transfer matrix eigenvalues, 
we derive the connections to the c—1 CFT together with the expressions for g. Since we find the FSS form for the 
general scaling levels, the TPF's for the periodic and the general twisted boundary conditions are derived. In Sec. 0, 
we summarize our results and discuss relations to other works. Also, possible physical relevances to the ECS are 
discussed. In App. |b| we prove an identity J — 1 where J appears in the relations between the finite-size correction 
terms and susceptibilities. In App. ^ the modular transformation properties of the TPF for the ASSV model with 
general boundary conditions are discussed. 

II. THE ASYMMETRIC SIX- VERTEX MODEL 

On the square lattice of N columns and of M rows, the six-vertex model configurations are obtained by covering 
the bonds of the square lattice with arrows which satisfy the ice rule: At each vertex there are two arrows in and 
two arrows out. The six vertex configurations satisfying the ice rule are shown in Fig. 1. Following the notation of 
Ref . [fol , the vertex energies are assigned as 



ei — —^ — h — v, £2 = —^ + h + v 

63 = ^ - h + v, £4 = I + (9) 
£5 = -e, £6 = -e . 

We write the energy in units of fc^T and denote the vertex weights as Wi = exp(— £i). h (v) is the horizontal (vertical) 
electric field which is conjugate to the horizontal (vertical) polarization 1 — 2pi (1 — 2p£)) where (pd) is the mean 
density of left (down) arrows. For later uses, we define parameters 

^ _ W1W2 + W3W4 - _ + - e^' ^^^^ 

2^WiW2W3W4 2 

and 

~ _ W1W2 + ■W3W4 - _ A ^^^^ 
W1W3 + cosh(2/i) 

The partition function Z is written as 

^ = Tr , (12) 

where T is the row-to- row transfer matrix whose definition can be found in Ref. |l0|| . Let Q denote the number of 
down arrows per row and let 
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The thermal average of 9 is pu. Since Q must be the same in the two adjacent rows of vertical bonds due to the ice 
rule, T can be considered in each subspace of fixed number of Q, which will be called the Q sector, separately. From 
the Bethe ansatz method m , the eigenvalue of T in the Q sector is given by 



A = A, 



(14a) 



where 



W2 



Wii — WiZj 
W4 — WiZj 



(14b) 
(14c) 



where the fugacities {zj} are given by the solutions of the Bethe ansatz equation 



(15) 



Note that the five-vertex model is achieved in the limit h — s- ±00 keeping A fixed. Introducing a variable p such that 
e'P = ze^'' and the phase function 



1 ^ 

^°,(p)-p + 2/iz--^eV 



(16) 



where 



1 + ^i{j>+v') _ 2Ae*P 
1 + e»(p+p') - 2Ae*P' ' 



(17) 



the Bethe ansatz equation takes the simple form 



2^ 

N 



(j = i,...,Q) 



(18) 



where Ij are half-integers (integers) for Q even (odd). If we define the energy and Eq through the relation 



Afl,L = exp 



v[N - 2Q) - Ej^-^ 



(19) 



they are given as 



(20) 



where 



^0 , ^ , 2A - e±* - e±'P 



(21) 



and {pj} are the distinct solutions of Eq. (|l§|). In Eqs. ( pO[ ) and (|2l|), the upper (lower) sign corresponds to the R 
(L) case. 

The bulk free energy is obtained from are ground-state energy. Using the ground-state energy, the free energy e as 
a function of h and q is given as 



e(/i, g) = hm min < — — , — — 



(22) 



4 



where the hmit is taken with fixed q. The free energy / as a function of h and v is given by the Legendre transformation 
of e, i.e., 



f{h,v)= mm{eih,q)^il-2q)v}. (23) 

0<<3<1 

The mean down-arrow density po is the value of q which minimizes the expression above while the left-arrow density 
PL is given by p^ — (1 + de/dh) jl. Of the four model parameters, only A and h enter the Bethe ansatz equation 
determining {pi} while $5? l depends only on A and 8. The vertical field v simply adjusts the mean value of q. Note 
that the horizontal field h plays a role of an additive constant to the phase function in Eq. (p^. This fact enables 
one to follow the same line of analysis as in The asymmetric XXZ Hamiltonian is diagonalized with the same 
Bethe ansatz method but with a different energy function . In next section, we will show that the FSS properties 
do not depend on the form of the energy function. So the asymmetric XXZ chain shares the same critical property 
with the ASSV model. 

Different choices of the set {/j} in Eq. (|l8|) lead to different eigenstates. It is well established that the ground-state 
energy is obtained if /j = — (Q + l)/2 + j. In analogy with the free fermion theory, we will say that a position j is 
occupied by a particle (hole) if = — (Q + l)/2 -I- j is included (not included) in the set {/j }. Then the solution is 
classified by the particle-hole configurations. The ground state corresponds to the Fermi sea as shown in Fig. 2(a). An 
important class of levels is characterized by shifting the Fermi sea by m, i.e., choosing Ij — <j)j + m. We call this the 
TO-shifted levels and show in Fig. 2(b) an example. Creating particles and holes at either ends of the m-shifted states 
generates the whole class of excited states which scale as 1/N in the critical phase. General form of these excitations 
is obtained by creating Up particles at positions j = Q + m + pk, {k=l,2,. . ., n^) with 1 < pi < P2 < ■ ■ ■ < Pn , 
holes at j — Q + m + I — hk, (fc=l,2,. . ., rih) with 1 < /ii < /12 < • • ■ < hn^, rip particles at j ~ 1 + m ~ pk, {k—1,2,. . ., 
fip) with 1 < pi < p2 < ■ ■ ■ < Piip , and fih holes at j = m + hk, {k—1,2,. . ., fih) with 1 < hi < h2 < ■ ■ ■ < hn^. 
Without loss of generality we can set Up = Uh and fip = nt- A particle- hole configuration is collectively denoted by V . 
And we will denote such an excited level by {Q,m,V). In Figs. 2(c) and (d), we give some examples of particle-hole 
configurations. 

For each range of A, there exists a variable transformation p = p{a) with which 0°(p(q;),p(/3)) defined in Eq. ( [1^ ) 
depends only on {a — /3) [p| pX|j2^ . The resulting function is denoted by 6, i.e., 6(0; — /3) = 9°(p(q;),p(/3)). Zt^(a) 
and ^R^h{oi) are defined similarly. We use superscript for functions of p and no superscript for functions of a. In 
the thermodynamic limit N ^ 00 with fixed g, the Bethe ansatz equation becomes an integral equation for the phase 
function under the assumption that the solutions of Eq. (|l^) lie densely on a smooth curve C in the complex-a plane 
with end points A and B. For the ground state, the contour is symmetric with respect to the imaginary axis so that 
A = —B*. Then Zoo = limjv— ►oo Zn should satisfy the following integral equation 

Zoo{a)=p{a) + 2hi-^ ( e{a~ P)Z'^{P)dp . (24) 

The solution of Eq. (U) depends on A and B which are determined as a function of q and h from the generalized 
normalization condition 

Zoo{B) = nq (25) 

or, equivalently, 

Zoo{A)^-TTq. (26) 
Using the solution of Eq. (|2^) and relation (|25|), the bulk free energy e in Eq. (|2^) is given by 

e(/i,g)=min|T/iT^-^^ $i^,L(a)Z^(a)da| (27) 

where the upper (lower) sign corresponds to the R(L) case. 

Before closing this section, we discuss the effect of the twisted boundary condition. The boundary condition do not 
affect the bulk properties but change the operator content and hence the TPF. The twisted boundary condition in the 
context of the XXZ chain is to impose the condition u^^i — exp(±27ri/)cr^ where af are the Pauli spin operators. 
In the vertex model language, it is equivalent to assign an extra vertex weight exp(±7i'iZ) to the horizontal arrow in 
the first column. More generally, one can impose the twisted boundary conditions (Z, I') by introducing the seams in 
the first column and the first row where an extra weight e~*'^' (e''^') is assigned to each right (left) arrow in the first 
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column and similarly an extra weight e"^' (e~"' ) is assigned to each up (down) arrow in the first row. Note that 
the effect of the vertical field is to give each up (down) arrow on the lattice an extra weight e" (e""). However, since 
the number of up arrows are conserved from row to row, one obtains the same effect if a vertical field of strength Mv 
is applied only to the first row of vertical bonds. Conversely, having a seam with extra weights e^*'^' is equivalent 
to assigning extra weights e*"' to all vertical bonds. A similar observation holds for the horizontal field also. 
Therefore, one then sees that the ASSV model with the twisted boundary conditions (Z, I') is equivalent to the ASSV 
model with periodic boundary conditions and with fields h = h — ttU/N and v = v + ttU' /M. In the XXZ chain, such 
relations are achieved by an unitary transformation as discussed recently in . A similar symmetry operation exists 
in the six- vertex model too pstl. 



III. SUMMATION FORMULA 



In this section, we derive the summation formula which converts the sum of the type 

Q 

N 



where /" is an arbitrary function which does not depend explicitly on h and q and {pj } is the solution of Eq. ( p^ ) 
for the {Q,m,V) level, into a series in 1/N. We obtain the expressions for the expansion coefhcients and derive the 
useful relations between the correction terms and partial derivatives of the bulk contribution of the sum. The phase 
function in Eq. (FgI) and the energy eigenvalue in Eq. (|2^) involve such kind of sums. For example. 



E^'^ = -N 



R,L\ 



With the change of variable p = p{a) explained in Sec. H, the sum becomes 



where /(a) = f^{p{a)) is not to be confused with the free energy f(h,v), aj are given by 



aj = 



1 f^TTlj 



N 



(28) 



(29) 



(30) 



with {Ij} corresponding to the level (Q, m, V), and is the inverse function of Zj^. Using Eq. ( |3C| ) in Eq. ( |29| ) and 
making explicitly reference to the locations of the particles, one can re-write Eq. (^) as 

-. ni 
k=l 

+ aFX! [/(^Ar^('/'l+™-pJ) - fi.Z^^i.^m+hk))\ 



k=l 



where 



Q + 1 



(31) 



(32) 



The first sum is over the Fermi sea, the second accounts for the shift of the Fermi sea, and the third (fourth) accounts 
for the particle-hole configurations at the right (left) end of the Fermi sea. The general strategy here is to regard 
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Z^^ as known first and determine it self-consistently later. To proceed, we make the crucial assumption that Zjj^{(j)) 
has the first derivative at (/) = ztTrg. We will see later that the scaling of the energy gaps, or the mass gaps, 
throughout the critical phase follows from this assumption. Conversely, we assume here that the critical phase is 
characterized by the fact that (iTrg) exists with possible exceptions at some special points such as at (/i=0, 

(7=1/2). Note that this assumption fails on the stochastic line A = 1 where the energy gap does not scale as 
but as l/Af3/2 [|. Now, applying the Euler-Maclaulin sum formula to the first sum of Eq. (plj) and using the Taylor 
expansion of f{Z^'^{(j))) at (j) — ±7rq for the rest of the sums, one then obtains to 0{1/N^), 



where 



1 

2^ 



7rq 



f{z^\^)) d<j>+-[f{z^'m] 



nq 
— TTq 



2. 




1 


iV2 




^ 24 


2. 


( ir? 


1 


iV2 


\~ 


^ 24 



M]r {Z-^^TTq)) Z-/ {nq) 



(33) 



fc=i 

rip 



(34) 



fc=i 



This is not a correct 1/iV expansion yet since Zj^'^{4>) also has TV dependence. We next assume that Zi^{a) possesses 
the 1 /N expansion of the form 



ZN{a) = Zoa{a) + 



hja) ^ b2{a) 
7V2 



N 



1 

iV2 



(35) 



where Z^o, and &2 should be determined later. Inverting Eq. 
dependence of Z^^{(j)) can be shown to be 



hia) 



Z'^ia) ZUa) \2Z'^ia) 



^) using the change of variable (j) — Z^ (a) , the N 
h\{a) 



1 



(36) 



where ' denotes the derivative with respect to a. Inserting Eq. ( p6| ) into Eq. (|3^), expanding all quantities to 0{1/N'^) 
again, and changing the integration variable from to a, one finally obtains the summation formula in the form 



S2_ 

iV2 



1 

N2 



where 







So = 


2n j 








Si = 


2^ J 








1 


52 = 






^2^ 



f{a)Z'^{a)da 



f{a) [2TTm - bi{a)] da 
— r f{a)b2{a)da + ^^^^^ 



(37a) 

(37b) 
(37c) 



^ 1- 



bi{B) 
27rm 



-24+^ 



2^r(A) 



ni 



1 - 



2Tim 



1 

24 



-AT 



(37d) 



Here A and B = —A* are determined from Eq. (|2^) or Eq. (^6|) as functions of h and q. 

In deriving Eq. (37), we assumed the scaling form of Zjv in Eq. (^5|). The self-consistent equations for Zqo, &i, and 
&2 are obtained by applying the summation formula to Eq. (16) with /(/3) = ~Q{a — (3). Equating corresponding 
orders in both sides of the resulting equation, we obtain Eq. (24) for Zoo (a) and 
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T o [bi{a)] = m[e(a -A)- e{a - B)] 
To [62(a)] ' 



ZLiB) 
2tt 



JTT / _ bi{B) 
2 \ 27rm 



-24+^ 



2 V 27rm 



24 



-AT 



K{a - B) 
K{a- A) , 



where K('-/) = dQ{'j)/d'y and T is a linear operator defined by 



r o [G{a)] = G{a) + — / K{a - (3)G{(3) d(3 
27r J A 

for any function G{a). 

The solutions of Eq. ( |3^ ) are written in terms of the function F{a, /i, A, B) which is defined by 

To[F{a,tJL,A,B)]=-]^Q{a- ^l) . 

Because Q is odd and A = —B*, it satisfies the relation 

Fia,^i,A,B)^^Fi-a*,-^I*,A,B) . 
Using the linearity of Eq. ( ^8| ) , 61 (a) can be written as 

61(a) = 27rm(l - D{a)) 

where 

D{a) = 1 - - [F{a, B, A, B) - F{a, A, A, B)] , 

TT 



(38a) 



(38b) 



(39) 



the dressed charge function Equation (41) implies that D{a) = D{—a*) and so D{A) = D{B). Usin; 
is written as 



(40) 

(41) 
(42) 

(43) 
this, 62(a) 



62(a) 



An f DqItl^ 1 



+ - 



ZUB) V 2 24 
47r ( Dqtu'^ 1 



Af] F {a,B,A,B) 



o 

where F (a, /x, A, B) = ~dF{a, fi, A, B) and 

Do = D{A) = D{B) 



2 -+JJ-] F{a,A,A,B) 



Inserting these into Eqs. (B7c) and (37d), we can put Si and ^2 in the form 



where 



Si = m f' {a)D{a)da , 
J A 

f Dh-n^ , 1 \ f Dlm^ 1 , 



/'(a) F (a, B, A, B)da + 7r/'(B) 



f'ia)Fia,A,A,B)da + nf'{A) 



(44) 

(45) 

(46a) 
(46b) 

(47a) 
(47b) 
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We have used the notation and C for the quantities on the right-hand side of Eqs. (47a) and (47t), respectively, 
anticipating identification of them as the anisotropy factor. The fact that C is the complex conjugate of C is not 
transparent in this form but will turn out to be the case as will be seen later. 

The resulting expressions for Si and 52 seem to be rather complicated. But, the following manipulations show that 
they are related to the partial derivatives of Sq with respective to q and h. To take necessary derivatives, one need 
to consider the variations of A , B, and Zoo (a), which is denoted as SA, SB, and dZoc{a), respectively, upon the 
variations of h and q, denoted by dh and Sq, respectively. We show in App. A that they are given by 



where 



From Eq. one can see t hat D 2(a) = —D2{— 

variation of Sq given by Eq. ( t37b| ) is 



Z'^{A) SA^ -{TT + D2{A))Sq-2iDoSh , 
Z'^{B)SB ^ {n- D2{B))Sq^2tDoSh , 
6Z^{a) = D2{a) 5q + 2iD{a) Sh , 



D2{a) = F{a, A, A, B) + F(a, B, A, B) . 

which implies that D2{A) 



(48a) 
(48b) 
(48c) 



(49) 



SSo 



-(fiA) + f{B))Sq- 



1 

2^ 



f'{a) SZoo{a) da, 



-D2{B). On the other hand, the 



(50) 



where a partial integration and Eqs. (Eq) and (E6[) are used. Combining Eqs. (pO[) and (48c), one obtains 



— / f{a)D{a)da 



dSo 
dh 

dSo 



1 



{I{A) + f{B)) 



which, compared with Eq. ( 46a ), gives the first relation 

Si = in 



1 

2^ 



dSo 



f'{a)D2{a) da 



dh 



(51) 
(52) 

(53) 



It means that the first-order correction term is proportional to the partial derivative of the bulk term with respect to 
h. Next, the variations of dSo/dh and dSo/dq are shown in App. |^ to take the form 



dSo 
dh 

dSo 



JiC + OSq- 



2iDl 



(C 



05h , 



2Dl 



where C,, C, are given in Eqs. (47a) and ( 47b ), and 

J = Z?o(l - D2{B)/7r) = Do{l + D2iA)/n) 



(54) 
(55) 

(56) 



J is a constant depending on h, q, and A through A, B, and 8. But we find that the value of J is equal to 1 for any 
values of A, and B only if the function is odd and A = —B*. We give the detail proof in App. Using the fact 
that J = 1, we then have from Eq. (|5j) the expression for C and C in terms of the partial derivatives of Sq as 



and 



c 



1 


" d 


fdSo\ 


ni d'^So 


1 


" d 


2 


dq 


\dhj 


^ 2Dl dh^ 


~ 2 


dh 


1 




rdSo\ 


ni d^Sa 


1 


d 


2 


dq 


ydh) 


2Dl dh'^ 


^ 2 


dh 



2iDl d^Sp 
TT dq^ 



dSo\ , 2^£)g d^Sp 
TT dq^ 



(57) 



(58) 



Since iSo is a real function of h and o, it is now obvious that C is indeed the complex conjugate of C as claimed before. 
From the second equalities in Eqs. (|57|) and (|8|), we find the identity 



9 



oq-' 



\2Dl 



(59) 



This identity comes from the fact that J 
phase of the ASSV model. Equations (|5^ 
the derivatives of Sq. An analogous relation between the Fermi velocity and magnetic susceptibility has been obtained 
by Bogoliubov et al. p4| for the symmetric XXZ chain but our result is more general. 



1 and is very important in identifying the universality class of the critical 
(p7|), and (|58|) together with Eq. (46b) relate the finite size corrections to 



IV. FINITE-SIZE CORRECTIONS OF THE TRANSFER-MATRIX SPECTRA 



In this section, we investigate the finite-size corrections to the transfer matrix spectra of the ASSV model using 
the summation formula obtained in Sec. [II. For simplicity we only mention the case where Eq < Eq which we call 
the R case. Most of results are the same for both R and L cases. The energy for the level (Q,m,V) is obtained by 
applying the summation formula to Eq. (p8|): 

EQ,^^-p = -(h+^-+So[^i,]^N-S,[^i,]-^^ + o(^^^ , (60) 

where the leading order term contributes to the bulk free energy, i.e., Sq — —e — /i — (5/2 from Eq. (p^). The correction 
term Si is obtained from Eq. ( |5^ ) 

Si = — 7rim(e/i + 1) (61) 



and ^2 is given by Eq. (46b) with 



^ = -7;(ih,q- ^^eh.,h = -^eq.h+"—eq.q (62) 

'lUn Z TT 



from Eq. (57). We use the short-hand notations that the subscripts of e denote partial derivatives. Note that 
e/i = 2pL — 1 where pL is the left-arrow density (For the L case, {eh + 1) in Eq. ( |6l| ) is replaced by {en — 1)). 
The partition function is 

Z ^ ^ (,-MEQ,„,,v+vMiN-2Q) ^gg^ 
Q,m,r 

with EQ^ra,!^ given in Eq. (|60|). To sum over the sectors Q, we take advantage of the fact that the summand in Eq. ( |6^ ) 
is peaked around the value of Q near Q = Npo- Inserting q = po + {Q — Q)/N to Eq. (|60| ) with a term —v{l — 2q)N 
added and expanding to order 1/A^^ assuming Q — Q is of 0(1), one gets the scaling form of the energy as 

EQ,„^.,v - v{l - 2q)N = /(/i, v)N + 2TnmpL - -jf- {m{Q - Q) + - 77) 



DW + - + + ^ - I^) + « (^) (64) 



where f{h,v) is given in Eq. (p3|), C' (C") is the real (imaginary) part of C, given in Eq. (|62|), and the relation 
dpL/dq — eh.q/2 has been used. 

On the other hand, for models whose TPF is given by Eq. (^, the 0(1/A^) part of the energy is expected to behave 

as 



27rC" (m^ gn^ —. l\ 2niC , ,r irr\ 



Comparing Eq. ( |6^ ) with this expression and identifying the indices m and Q — Q in Eq. (^) as the spin-wave and 
vortex indices, respectively, as in the SSV model, one sees that the critical phase of the ASSV model is indeed in the 
Gaussian- model universality class with c = 1. The Gaussian coupling constant g can be read off from the coefficient 
of in Eq. (64) as 

^ = 2^ 
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where we recall that Dq is defined in Eqs. ( |45| ) and (^). However, from the coefficient of (Q — Q)^ in Eq. ([6^), we 
also have the relation 



Therefore, to identify the spectra Eq. (|6^ ) with those of the Gaussian model, we require that the two expressions of 
g, Eqs. (|65|) and (^6|), give the identical result. Using Eqs. (|62|) and (|65|), the condition is then 



^ = {^gf (67) 



which is guaranteed by the identity (g9|) . That the coupling constant g is given by Eq. ( |65[ ) is first derived by Izergin 
and Korepin for the symmetric XXZ chain in a magnetic field but it also holds for the ASSV model. Using the 
relations 

fh,hfv,v fh^vfv.h 

S-h,h = 7 — 

Jv,v 

4 

Jv,v 

in Eq. (|67|), we obtain an alternative expression given in Eq. (|7|). 

Before proceeding further, we should ask the validity of the summation formula. If the summation formula is valid 
with non-trivial solution for i-e., &2(q^) 7^ 0, then transfer matrix spectra follow the scaling form in Eq. which 
holds only in critical phase. We have assumed the expansion forms of Zjq and Z^^ in Eqs. (^5|) and ( |36| ) in deriving 
the summation formula. The assumption is valid only when the coefficients are finite, i.e., Z'^ should not be 0. 
Otherwise, the summation formula does not hold and the criticality is not guaranteed. So it is the sufficient condition 
for the criticality that Z'^{a) and 62(0;) are not zero. It is compatible with earlier studies on the phase diagram of the 
asymmetric six- vertex model. First, for A > 1 Z'^{A) is always for any values of q when A = 1 ^J. At this point, 
the ASSV model describes stochastic growth models and the finite-size corrections show the Kardar-Parisi-Zhang 
scaling [Q,^. Another case occurs when A < — 1. In this region, there exists the anti- ferroelectric phase when q = 1/2. 
The phase boundary is given in ||l^ and Z'^{A) at the phase boundary is also 0. The finite-size corrections on this 
phase boundary has recently been studied by Albertini et al. p^ . There are also ordered ferroelectric phase with 
q = in which case A = B from Eqs. (|25| ) and (^). In this case 62(0;) is identically zero and the system is also out of 
criticality. The q — and 1/2 ordered phases are separated from the critical phase by the Pokrovsky-Talapov (FT) 
transition 

Inside the phase boundaries stated above, t he s ystem is critical and the finite-size corrections to the transfer matrix 
spectra are obtained from Eqs. (|60|), (|6l|), ([46b| ), and (|6^). Actual values of g through the critical phase can be 
mapped out without difficulty. In Figs. 3, we present the constant-g lines in the (A, q) plane for tanh(2/i) = 0, 0.5, 
and 1.0. As seen from Eq. (|65|), g is obtained from the function D. Instead of solving Eq. ( ^0| ) for F and using 
Eq. (^), we study the integral equation for D{a) 

To[Dia)] = l, (68) 

which is obtained by combining Eqs. (^3|) and (^o|). It can be solved analytically only for special cases. For q — 0, 
there is a trivial solution Do{a) = 1. So we have g — 1/2 at (7 = and at g = 1 by symmetry. For A(= — cosh A) < —1, 
the anti-ferroelectrically ordered phase with g = 1/2 appears for A < Ac with 



where he is given by [nol 



he = In cosh A V ^e~2"^ tanh(An) . (70) 

n=l 



The value of Ac decrease from —1 for tanh(2/i) = to —4 for tanh(2/i) = ±1. At the phase boundary, the value of B 
is TT-f iA [0. Then Eq. ( |68| ) is solved by Fourier series method to yield the solution D(a) — 1/2. So we have 5 = 2 at 
the anti- ferroelectric phase boundary. For A > 1, the ASSV model is critical only for A < 1. As A approaches 1 from 



11 



below, g decreases to 0. For the SSV model, the model is critical for —1 < A(= — COS7) < 1 and the value of g is 
given by the simple formula g = 1 — ^ /tt jl^. In this case, B ^ 00 and naive application of Eq. (|68| ) is invalid due to 
the fact that Z'^{±oo) = 0. However, a careful treatment using the Wiener-Hopf method can reproduce the correct 
result starting from Eq. (|8|) For other cases, the values of g are obtained by solving the integral equation 

( |68| ) numerically. Note that there is a discontinuity in g from 1 to 2 at (A, q) = (—1, 1/2) for tanh(2/i) = 0. It is the 
Kosteritz-Thouless (KT) transition point where the free energy has the essential singularity . There is a crossover 
from the KT transition to the PT transition for tanh(2/i) 7^ 0. 

We next derive the TPF Z. Using the scaling form of Eq. (p4) and the identity Eq. (|67|), it can be put in the form 



2TTiT A 



^m,Q-Q 



2mf (A ^ t^+AT- — 



(71) 



where the modular ratio r is given as r = ^C- The sum over V produces the factor 



1 



'?(q)'7{q) 



where the nome q — e^'^" and 7^(q) is the Dcdckind eta function given in Eq. (|). Now we introduce the two mismatch 
parameters 



= {plM} and 13 = {poN} = {Q} 



(72) 



where {x\ denotes the fractional part of x\ i.e., a and /? are the mean number of left arrows per column and down 
arrows per row, respectively, modulo 1. We next change the index Q in Eq. ( |7l] ) to n by the relation Q — Q = n — f3. 
Then the TPF for the critical phase of the ASSV model is given in the form of the modified Coulombic partition 
function as given in Eq. (^). The mismatches come from the incommensuration of arrow densities pL and pjj with 
system size M and N, respectively. 

We also consider the effect of the twisted boundary condition on the TPF. Since the twisted boundary conditions 
{I, I') corresponds to the ASSV model with the periodic boundary conditions but with the modified fields h — nil/N 
and V + nil' /M, these 0{1/N) changes of h and v modifies the TPF. Thus, one replaces h and v by h — nil/N 
and V + nil' /M, respectively, in Eq. (p4), and makes Taylor expansion to the necessary order in 1/A'^ to obtain the 
boundary condition effect on the finite-size corrections. After a straightforward calculation the resulting TPF is found 
to be generalized to 



ZipL,pD; I, I'; M, N\t) = e^'*^'(2''--i)-"^''(2p--i)Zgc({Afp4, {Npoh I, l'\r) 
where Zgc is the generalized Coulombic partition function defined as 



Zgc (a, /?;/,/' |t) 



h(q)P 



E 



^-27Tima-27Til' {n-j3) ^Am-l,n~f3 ^^m-l,n-l3 



(73) 



(74) 



The toroidal partition functions in Eqs. (|§) and ( [73|) satisfy the necessary modular covariance as they should be. 
These are discussed in App. ^. 



V. DISCUSSIONS 



In this paper, we utilized a method of calculating the finite-size corrections applicable to the critical phase of the 
ASSV model and investigated the ESS of the transfer-matrix spectra of the model for arbitrary sector Q — qN. Since 
we work for general the string solutions appearing only at g = 1/2 do not complicate the analysis and low- lying 
excitations are easily classified. For any low-lying levels, the finite-size corrections for arbitrary sum over functions of 
rapidities are related to partial derivatives of the bulk contributions; Eq. ( p3[ ) for the first order term and Eqs (46h), 
{pT(}, and (|5^ ) for the second order term. From this, finite-size scaling amplitudes of the ASSV model are related to 
thermodynamic quantities. 

The ASSV model is conveniently parameterized by the interaction parameter A (Eq. ([ll|)), the anisotropy S 
(Eq. (^)), the horizontal field h and finally vertical field v or alternatively, the down arrow density q. For fixed 
h, S, and q{^ 1/2), the model is critical for —00 < A < 1. When q = 1/2, however, the region defined by A < Ac 
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with Ac given by Eq. (|6S|) corresponds to the anti-ferroelectrically ordered phase. We showed that the critical phase 
is in the Gaussian-model universality class with c = 1 and the Gaussian coupling constant g as given by Eq. (|6^). 
The value of g ranges from at the stochastic limit A = 1 to 1/2 at the PT transition lines bordering the g = or 
q = 1 phases and to 2 at the PT transition lines bordering the q = 1/2 anti-ferroelectrically ordered phase. 

We also constructed the TPF in Eq. for periodic boundary conditions and in Eq. (|7^) for more general twisted 
boundary conditions. Due to the incommcnsuration of arrow densities, it takes an additional phase factor e"^'^""" 
and a constant shift (3 in the vortex-excitation index n. The phase factor comes from the 0(1) imaginary term of the 
energy spectra Eq. {pW which can be understood as follows. The continuous lines connecting the left pointing arrows 
may be pictured as domain walls or steps running across the lattice. Since there are p^M of them, the mean distance 
between them is 1/pl lattice units. Thus the periodicity of the lattice in the time direction is enlarged by the same 
factor and the transfer matrix T is the p^-th power of the time translation operator for one unit distance. Therefore, 
a factor exp(27ri/OL) may appear in the spectrum. As shown in App. |c[ the mismatch parameter in the phase factor of 
Eq. (H) becomes that which shifts the vortex-excitation index upon the modular transformation which simply rotates 
the lattice by 90° exchanging the role of M and TV. 

Our result for the TPF is consistent with all other previous results. For the five-vertex model in the noninteracting 
case, the TPF on a deformed square lattice is given as 

^5-v = ^ J2 e--'^«^'"q^"-"-«oq^— «o , (75) 

where Qo = {1 — Pd)N and Qi = (2pi, — pn — 1)M. The TFP's on different geometry can be transformed to each 
other through the modular transformation [See App. |c|. One can find that Eq. ( [75|) is obtained from Eq. (^) through 
the modular transformation r t + M / {2N) with M/(2N) — 1. In the SSV model, pL—p£i~\/2 and hence a {(3) is 
if M (N) is even and 1/2 if odd. If the twisted boundary condition (1,0) is applied in the SSV model, the operator 
content of this model obtained by other methods is summarized by the TPF 

^=u;7^ E q^"— "q^"— " (76) 

where /i(i/) is for M(N) even and 1/2 for odd. It is in agreement with Eq. ( [73| ) when one uses Pl = Pd = h and 
/' = in Eq. (|7|). 

The asymmetric XXZ chain is obtained if one takes an extreme- anisotropic limit in which the vertex-weights 
approaches the limit = w/^ = Q and wi = W2 = = wq = 1. With an appropriate parameterization, T can be 
put in the form 

T(u) =exp(zP-7.Hxxz+0(^')) (77) 

where u is the so-called spectral parameter as a function of which T(u) form a commuting family, P is the shift 
operator, and Hxxz is the asymmetric XXZ chain Hamiltonian. Therefore, T and Hxxz share a same set of spectra. 
However, the TPF is not the same as Tr exp(— uMHxxz)- This is why complete information on the operator content 
for the XXZ chains is not sufficient to construct the TPF of the lattice model. 

When the ASSV model is considered as a model for the ECS of e.g. fee (110) surface, the free energy /(/i, v) itself 
is the height of the surface from the base (110) plane with appropriate identification of coordinates. In particular, the 
surface curvature n defined in e.g. 1^ is related to the Hessian of / by 

n = ^H^'^ (78) 

where d is the distance between the crystal planes. Our result Eq. (^) then relates k to g by 

2cP 1 

K — — ; 

ttkbT g 



(79) 



exactly. When the interactions are intrinsically antiferroelectric, the q = 1/2 phase at low temperatures corresponds 
to the flat (110) facet. As the temperature is raised the facet area becomes smaller and finally disappears at the 
roughening temperature Tr. The roughening transition corresponds to the KT point where g = 1. Thus there occurs 
a universal jump in the curvature of magnitude 2cP /(irkBTfi). This universal jump has been anticipated from the 
Coulomb gas picture of the solid-on-solid models Q and also from the solution of the six- vertex model near h=v~0 ^ 
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and has been measured experimentally ||3^]. Our result goes beyond this. Equation ( [79| ) being exact, it applies to 
all points of the curved portion of the crystal. In particular, near the PT transitions where the curved surface joins 
the facet smoothly with the exponent ^ = 3/2 , one expects the universal jump of the curvature with g — 2 or 
g = 1/2 depending on which portion of the PT line is appropriate. It is interesting to speculate whether the relation 
Eq. ^79) also holds for other solid-on-solid type models and for real samples. If it is valid for real systems, it would 
provide a mean to measure the critical exponent directly from surface curvatures. 
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APPENDIX A: VARIATIONS OF THE PARTIAL DERIVATIVES OF 5( 



In this appendix, we consider the variations of Sq upon the variations of h and q. As h and q vary, the end points 
A and B also vary satisfying 



TrSq = SZ^{B) + Z'^{B)SB 



(Al) 



and 



-7:Sq^6Z^{A) + Z'^{A)6A . (A2) 

They are obtained from the variations of Eqs. ( psj ) and (p6|). Since this involves dZoc{a), we take the variation of 
Eq. (^) after performing the partial integration inside the integral operator. This gives a simpler result 



SZooia) 2iDia) 6h + Da (a) Sq 



(A3) 



with D and D2 defined by Eqs. ( ^ ) and (^). In taking the variations in this appendix, it is convenient to use the 
integral representations of D and D2 rather than the definitions in Eqs. ( ^ ) and (|4^). The integral representation 
for D is given in Eq. (^8|) and the other for D2 are easily obtained by combining Eqs. (|4|) and (^0|): 



To[D2ia)]^--[e{a~A) + e{a-B)] 



where the operator T is defined in Sec. ttl|. Inserting Eq. (A3) into Eqs. ( [Al[ ) and (A2) gives Eq. (| 
Next, consider the variation of Eqs. ( |51| ) and (^2|): 



(A4) 



dh 



TT 



f{a)5D{a)da + f{B)D„5B - f'{A)D„SA , 



fiA) ^ D2{A) 



5 A 



D2{B) 



SB 



(A5) 



(A6) 



The variations of D(a) and -02(0;) are obtained from Eqs. (|68| ) and (A4). After a straightforward calculation, one 
can find that 



SD{a) = — |f (a, S, A, B) SB- F (a, A, A, B) m| , 



5D2{a) = - 1 - 



D2{B) 



F {a,B,A,B)5B - 1 



D2iA) 



F {a,A,A,B)SA . 



(A7) 
(A8) 



Using Eqs. ( [A7| ) and ( [AS| ) in Eqs. ( [A5| ) and ( [Aq ) and using the definitions of ( and ( in Eq. ( |47[ ), one obtains 
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dSo 
dh 



{tCZUB)6B-zCZUA)6A) 



D2iA) 



(A9) 
(AlO) 



which together with Eq. (48) give Eqs. (|54|) and (55) 



APPENDIX B: THE PROOF OF J 



In this appendix we present some identities between the partial derivatives of F{a, fi, A, B) defined in Eq. 
prove the identity J = D{B){1 — D2{B)/'k) = 1, where 



and 



D[B) - 1 - -[F{B, B, A, B) - F{B, A, A, B)] 

TT 

D2{B) - F{B, B, A, B) + F{B, A, A, B) . 



(Bl) 
(B2) 



In general, J may be a function of a and b with A = —a + ih and B = a + ib. But, Eq. ( pO[ ) is invariant under the 
shift of all arguments by an imaginary amount, i.e., F{a, ^, A, B) = F{a + iw, fx + iu, A + iu, B + iu) for any real m. 
This means that D{B), D2{B), and thus J are functions of only a. When a — the function F(a, /i, A, B) becomes 
simplify —^<d{a — /i) and so J = 1. Thus, for the proof of J = 1, it suffices to show that the total derivative of J 
with respect to a is identically 0. 

The total derivative of J contains the partial derivatives of F{B, B, A, B) and F{B, A, A, B) with respect to a: 



da 
d_ 

da 



F{B, B, A, B) = F'- F -Fa + Fb , 
F{B, A, A, B) = F'+ F -Fa + Fb , 



(B3) 
(B4) 



where F' = dF/da, F= —dF/dn, Fa = dF/dA, and Fb = dF/dB. We use the convention that the arguments of 
the functions are (a, /i. A, B) when not shown explicitly. First, we derive several identities between partial derivatives 
of F, which will simplify Eqs. ( |B^ ) and (B4). Taking partial derivatives of Eq. (^0|), one obtains integral equations 
for F' , Fa, and Fb- After some manipulations, they take the form 



T oF' = --K{a - + — [F{B, /i. A, B)K{a - B) - F{A, /i. A, B)K{a - A)] 
2 27r 

ToFa^ ^F{A, A, B)K[a - A) , 
r o Fb = -^F{B, ^JL, A, B)K{a - B) , 



(B5) 
(B6) 
(B7) 



where T is defined in Eq. (p9|). Note that F and F can be written as 



F{a,iJi,A,B)^r-^o[--Q(a-ii)] , 
F (a, A, B) =. o i-^Kia - fi)] , 



where T is the inverse operator of T. Using the linearity of Eqs. (B5) 
obtained by acting T~^, which yields 



(B8) 
(B9) 



), and (B7), F' , Fa, and Fb can be 



F'(a, fi, A,B)^F {a, A,B)-- F{B, fi. A, B) F {a, B, A, B)- 

TT L 

F{A,^i,A,B) F {a,A,A,B) 
FAia,fi,A,B) = --F{A,^i,A,B) F {a,A,A,B), 

TT 

FB{a, fi. A, B) = -F{B, fi. A, B) F {a, B, A, B) . 

TT 



(BIO) 
(Bll) 
(B12) 
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Using these relations in Eq. (B3) and (B4), one can obtain 



-^F(B, B, A, B) = -F{A, B. A, B) F (B, A, A, B) 
da TT 



—F{B,A,A,B)^2F{B,A,A,B) 
da 



1 + -F{A, B, A, B) F {B, A, A, B) 



from which the derivatives of D{B) and D2{B) can be written as 

^1 P(B,A,A,B)ll + - [F{A, A, A, B) - F{A, B, A, B)] 
da IT [_ TT 

^5^1^ = 2p{B,A,A,B)ll + - [FiA, A, A, B) + F{A, B, A, B)] 
da TT 



(B13) 
(B14) 

(B15) 
(B16) 



As explained in Sec. F{A, A, A, B) = -F{B, B, A, B) and F{B, A, A, B) = -F{A, B, A, B) because the function 
G is odd and A = —B*. Using these properties in Eqs. (B15) and (B16), one can find that 



±D{B) = - °F{B,A,A,B)D{B) 
da TT 



d r D2{B] 
da \ TT 



-F{B,A,A,B) (l--D.2{B) 

TT \ TT 



(B17) 
(B18) 



Then the total derivative of J is obtained straightforwardly. 



dJ_ _ dD{B) / lMBy\ ^ jjfQxA. I 1 
da da \ / da 

= . 



D2{B) 



(B19) 



Equation. (B19) together with J = 1 at a = implies that J = 1 for any value of a. 



APPENDIX C: MODULAR TRANSFORMATION PROPERTIES OF Z 

In this Appendix, we examine the modular covariance of the TPF for the ASSV model in Eq. (|73|). It is convenient 
to work with the Gaussian model TPF with the shift boundary conditions fl^ : 



where 



-Zm,n(T) = Zoexp -^|n-rm| 



(CI) 



(C2) 



t" ?7(q)?y(q) 

is the Gaussian model TPF with the periodic boundary conditions. They have modular transformation properties as 

Z'm.niT) — ^■m,n+m{T + 1) — -2_„ ) (C3) 

or combining these 

' ar + 



(C4) 



where a, 6, c, and d are integers with ad — be = 1. 

The generalized Coulombic partition function defined in Eq. (Q) is written in terms of Z„i_„ by applying the Poisson 
sum formula: 



g{2TTm) 



(C5) 



m— — oo 



rn— — oo 
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with 

dxf{x)e'''- . (C6) 

-oo 

First, Eq. ( |7^ ) can be written as 

Zgc(a, /3; I, 1'\t) = ^ exp [-TTgT"{n - Pf - 2ml' {n - f3) - 2mal\ 



X ^ exp 



- — {m-lf + 2m{T'[n - [3) - a){m - I) 



(C7) 
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Then, applying the Poisson sum formula to the sum over m in Eq. ( |C7| ), one can express Zgc in terms of 2^n,m as 

Zgc(a,/3;Z,Z'|T)=Zo(r) ^ e-2-(™+")'-2-^''("-«Z„_^,„,+„(T) . (C8) 

m , 71 £ Z 

And the TPF in Eq. (|7|) becomes 

Z{pL.PD]l, I'] M, N\t) = exp [TnMli2pL - 1) - T:iNl'{2pD - 1)] Zo(t) 



X 

m , n G Z 



where a = {AIp^} and /? = {Npo} where {x} denotes the fractional part of x. 
Under the transformation t r + 1 , 



m , n G Z 

= Zo(r+l) ^ e-2"'("+"-/')-2"(''-0("-/3)z„„^,„+„_^(r + l) 



m , n G Z 

= i'gc(a -/?,/?; r - Z|t + 1) , (CIO) 

which means that 

Zgc{a,f3;l,l'\T+l) = Zgc{a + f3,f3;l,l' + 1\t) . (CU) 
So, the TPF Z transforms under the transformation r ^ t + 1 as 

ZipL,PD;l. I'; M, iV|T + 1) - Z{p'l,pd]1, l' + l;M + iV, A^jr) (C12) 

where 

, AfpL + Npp 
^ M + N ■ 

Similarly, Zgc transforms under the transformation r ^ — l/r as 

m.nGZ 

= Z^c{-P,a;l',-l\-i) , (CM) 

which means that 

Zgc{a,f3-l,l'\-L) ^ Zgc{f3,-a;-l'J\T) . (C15) 
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Prom this, one can also find that 

Z{pL,PD\ I, I'; M, N\-L) = Z{p'L,p'n; -I', I; N, M\t) (C16) 

where 

p'l=Pd, p'd = 1-Pl. (C17) 
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FIG. 1. Six vertex configurations. 
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FIG. 2. A few examples of particle-hole configuration for N — 24 and Q = 12. We denote the values of Ij by the angular 
coordinate of equi-spaced small circles on a unit circle. Closed circles denote the occupied positions, (a) depicts the the ground 
state, (b) the m — 2-shifted state, while (c) and (d) shows excited states from (a) and (b), respectively. 

FIG. 3. The constant-g lines are plotted in the (A,q) plane for tanh(2/i) = in (a), 0.5 in (b), and 1 in (c). These are 
obtained from the numerical solutions of Eq. (^8|). On A = and g = lines, g = 1/2. On A = 1 line g = 0. The value of g 
increases in steps of 0.1 from right to left and reaches 2 at g = 1/2 and A < Ac. The figures at the top show g as a function 
of A at g = 1/2 and the figures at left show g as & function of q at A — —3 in (a), —4 in (b), and —5 in (c). The rectangular 
symbols denote the phase-transition points Ac from the critical phase to the g = 1/2 ordered phase. 
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